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O Abstract 



. We present a general method for determining screened Coulomb parameters in molecular assemblies, in particular organic molec- 
C^l ular crystals. This allows us to calculate the interaction parameters used in a generalized Hubbard model description of correlated 
f^i organic materials. In such a model only the electrons in levels close to the Fermi level are included explicitly, while the effect of all 
D other electrons is included as a renormalization of the model parameters. For the Coulomb integrals this renormalization is mainly 
due to screening. For molecular materials we can split the screening into intra- and inter-molecular screening. Here we demonstrate 
£f) how the inter-molecular screening can be calculated by modeling the molecules by distributed point-polarizabilities and solving 
the resulting self-consistent electrostatic screening problem in real space. For the example of the quasi one-dimensional molecular 
i— h metal TTF-TCNQ we demonstrate that the method gives remarkably accurate results. 
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1. Introduction 

Low dimensional organic molecular compounds based on the 
BEDT-TTF, TTF and TCNQ molecules are narrow band sys- 
tems, where the Coulomb repulsion is large compared to the 
bandwidth. This leads to many-body effects which are man- 
ifest, e.g., in the spin susceptibility [1] or angular-resolved 
photoemission experiments [2, While theoretical works 

suggest the importance of Coulomb parameters in these mate- 
rials 01 [H, only rough theoretical estimates of these param- 
eters do exist 16D. The use of modern techniques, such as 
density-functional theory (DFT), gives accurate results of the 
parameters for individual molecules |[7U8t], improving previous 
quantum chemistry calculations lldll . However, there is a 
lack of accurate theoretical estimates of the Coulomb repulsion 
screened inside the crystal. These are needed for a more realis- 
tic description of organic molecular compounds. 

In the present paper we introduce a systematic method for 
calculating these Coulomb parameters in low-dimensional or- 
ganic molecular crystals. In Sec. [2] we construct the minimal 
model and give details of the general formalism. This is applied 
to TTF-TCNQ salts in Sec. [3] where we discuss the distributed- 
dipole approach. In Sec. [4] we give our results for Coulomb 
parameters. 



2. Model and Formalism 

The minimal model generally used to describe the low- 
energy electronic properties of low-dimensional organic sys- 
tems is the Hubbard model [11], where transfer integrals be- 
tween neighboring sites and on-site Coulomb interactions are 



taken into account. Such a simple Hubbard model is, however, 
usually not sufficient to describe strongly correlated organics, as 
long-range Coulomb repulsion energies can rarely be neglected 
10, Q • In me generalized model 
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we include hopping integrals between nearest-neighbors, f, on- 
site Coulomb energies, U, and longer-range Coulomb interac- 
tions, Vij. While transfer inte gral s can be obtained from ab 
initio DFT calculations J3,[l3, 14], Coulomb parameters are 
more difficult to obtain, as they include the screening correc- 
tion due to different high-energy processes inside the crystal, 
the main mechanism being the polarization of the surrounding, 
usually highly polarizable, molecules. By realizing that organic 
molecules preserve their identity inside the crystal, we may es- 
timate the screening energy, by describing the molecular inter- 
actions through classical electrostatics. 

In BEDT-TTF and TTF-TCNQ families of organic conduc- 
tors 11511 . flat molecules crystallize forming low-dimensional 
anisotropic structures. Considering such a lattice of organic 
molecules, we ask what happens when a charge is added to one 
molecule. The electric field of that charge will induce a dipole 
moment in the surrounding highly polarizable molecules. Solv- 
ing this electrostatic problem assuming fixed polarizabilities, 
we obtain the screening of the charged molecule in linear re- 
sponse. This gives the screening contribution of the lattice that 
is needed for renormalizing the bare Coulomb integral. 

Modelling each molecule as a polarizable point with charge 
q„ dipole moment p, and polarizability a, the energy of the 
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lattice in the above situation is given by 
"IP/I 2 
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+ ^ [C7i( P/ , pj) + U 2 (Pi, qj) + U 3 (g h qj)] , (2) 
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where E"' is the electric field and (£>(r) the potential associated 
with the external charge, while the first and fourth terms corre- 
spond to the resistance to polarization of each lattice point and 
to the internal energy of the molecules as a function of charge 
state, respectively. The terms U\, U2 and U3 are the dipole- 
dipole, dipole-monopole and monopole-monopole interactions 
II 1611 . For molecules without net charge (the electrons involved 
in charge transfer are explicitly included in the Hubbard model), 
we can simplify Eq. ©. Expressing it in the notation of Allen 
Intl . defining the dipolar field by \E dip ) = T\p), where T is the 
dipole-dipole interaction matrix: 



i + j, fi,v = x,y,z, (3) 



the energy of the system is then given by 
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-(pK 1 -r|p)-< P |E"'). 



(4) 



By applying the variational principle, SW = 0, we obtain the 
expression for the set of dipoles which minimize the energy of 
the whole lattice, 



l P > = (a- 1 -r)" 1 iE^>. 

Inserting Eq. (f5]l into Eq. © yields 

SW = —{E ext \ (a' 1 - |E"'). 



(5) 



(6) 



Therefore, placing two charges at a lattice point or two differ- 
ent lattice points, we are able to calculate the correction to the 
total energy due to the polarization of the rest of the molecules 
in the crystal. Defining the external field as the composition of 
two fields, corresponding to charges placed at points n and m in 
the lattice, E ext = E„ + E m , the energy of the system is given by 



anisotropy in the dipole interactions. Such instability influences 
the screening of a charge in the lattice. For small polarizabili- 
ties a the induced dipoles simply arrange along the field-lines 
of the central point charge. Increasing a we, however, observe 
quite unconventional dipole arrangements as shown in Fig. Q] 
This reflects the fact that for this lattice the matrix (or 1 - T) 
becomes singular at a critical polarizability a c , with the eigen- 
vector of the vanishing eigenvalue having non-zero momentum 
(in this case (n, n, 0)). Hence the ferroelectric instability hap- 
pens before the critical value acM obtained from the Clausius- 
Mossotti relation iv. 



3. Submolecular approach 

In a real molecular crystal the molecules typically form 
stacks, i.e., they are quite closely packed. Consequently, ap- 
proximating the polarizability of a molecule by just one point- 
polarizability is not a good approximation. It only works out- 
side a sphere containing the charge density of the molecule, 
which for typical 7r-bonded molecules will intersect the neigh- 
boring molecules, but it becomes accurate at large distances. To 
obtain a better description of the molecular response, we model 
it as point-polarizabilities distributed over the non-hydrogen 
atoms of the molecule B19I1 . The accuracy of this submolecular 
approach is tested for TTF-TCNQ salts. 

We consider a small system of only three or five molecules 
in a TTF chain. For such small systems we can perform con- 
strained density-functional calculations U3, |2QU, constraining 
the added charge to the central molecule. This gives us the 
screening due to the molecules in this small assembly. We 
compare to the results of the screening calculated using the 
distributed point-polarizability approach. Since this approach 
works least well in the near-field, this is a critical test of the 
method. Nevertheless, both approaches agree within 5%. If we 
included more distant molecules, for which the computational 
cost of constrained density-functional calculations would be 
prohibitive, the agreement would improve further, as the multi- 
pole approximation becomes exact in the far-field. Thus, using 
polarizability tensors obtained with DFT [Q] in our submolecu- 
lar approach essentially reproduces the results of full quantum 
mechanical treatment of the screening. 



SW = -<E m | (aT 1 - r) ' |E„) - <E„| (aT x - r) ' |E m ). (7) 4. Screening parameters in TTF-TCNQ 



The equation above gives the on-site Coulomb screening, SU, 
when n = m, and the inter-site, SV, SV, for n + m. Thus, 
defining \n - m\ = I, the general equation for the screening 
Coulomb parameters reads 



<5V' = -<E I (or 1 -r) -1 iE z > 



v l = V l + SV 1 , 



(8) 



and V l are 



where 6U corresponds to / = 0, 6V to I = 1, 
the Coulomb parameters for a single molecule/dimer, calcu- 
lated with DFT |0 H . 

To illustrate the method we apply it to a cubic lattice of po- 
larizable point dipoles, in which an antiferroelectric instabil- 
ity occurs above a critical polarizability a c IU7U18I1 . due to the 



To obtain the screened parameters for a TTF-TCNQ crystal, 
we perform calculations for crystal fragments (clusters) of in- 
creasing sizes to extrapolate to the infinite-size limit. To obtain 
the parameters for the extended Hubbard model ([TJ we consider 
the Coulomb screening between two charges at the same or at 
neighbor molecules placed in the center of the cluster. For sim- 
plicity, the clusters of molecules are constructed as spheres of 
radius R, where R is the distance between the doped molecule 
(or center of the pair of molecules) and the farthest neighbor 
in the cluster, hence Q = being the volume of the sphere 
1 2 1 ] . For large enough R we can linearly extrapolate the screen- 
ing energy 6U or SV, respectively, versus 1 /R to the thermody- 
namic limit (R — » 00). This is shown in Fig.|2l where clusters of 
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Figure 1 : Dipole arrangement in the x-v-plane of a simple-cubic lattice of point dipoles to a charge added at the center. As the polarizability increases from 
a = 0.5aciu (left), 0.75acM (middle) to a = Q.SacM (right), the dipoles start to deviate strongly from pointing along the field-lines of the external field. 
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Figure 2: Screening of on-site Coulomb interaction, SU, in a TTF molecule 
(left panel), and screening of inter-site Coulomb interaction, SV, for nearest- 
neighbor TTF molecules (right panel), versus fi~'' 3 , for clusters of increasing 
radius. Extrapolation to the thermodynamic limit is shown by the dotted line. 



SU 6V 6V 6V" 



TTF 
TCNQ 



2.7 
2.6 



1.9 
1.9 



1.25 
1.3 



0.8 
0.9 



Table 1: Screening Coulomb energies for TTF-TCNQ. SU is the on-site screen- 
ing and 6V, SV and SV" correspond to first, second and third nearest-neighbor 
screening energies between molecules in the same TTF or TCNQ stack. All 
energies are in eV. For comparison, the band-width is about 0.7 eV. 



up to 400 molecules are studied. Calculations for TTF-TCNQ 
salts up to the third nearest-neighbors are performed and results 
are given in Table \T\ The ratio of the radius of the sphere, R, 
to the distance between third nearest-neighbors, d, where we 
still can make the extrapolation to the thermodynamic limit is 
d/R = 0.43. However for longer range interactions, the cluster 
size should be increased towards this convergence ratio, reach- 
ing a linear dependence with the inverse of the volume, which 
allows the extrapolation to the bulk value. 



to other crystals like quasi-bidimensional salts based on the 
BEDT-TTF molecule. 
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5. Conclusions 



We have presented a general method for obtaining accu- 
rate values of Coulomb parameters in organic molecular com- 
pounds, where molecules interact weakly, and shown how it 
works for TTF-TCNQ. We are currently extending this method 
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